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ABSTRACT 


Starting  from  the  Painleve-Backlund  equations  obtained  from  Painleve 
analysis  of  the  Korteweg-de  Vries  equation,  closed  form  solitary  wave 
solutions  are  explicitly  constructed.  It  is  shown  that  repetitive  application 
of  the  Mobius  group  of  fractional  linear  transformations  does  not  lead  to  new 
solutions.  Various  connections  of  the  Painleve  method  with  Hirota's 
formalism,  the  Backlund  transformation  method,  the  Lax  approach  and  the 
Inverse  Scattering  Technique  are  discussed  in  detail, 
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CONSTRUCTION  OF  SOLITARY  WAVE  SOLUTIONS  OF  THE 
KORTEWEG-DE-VRIES  EQUATION  VIA  PAINLEVE  ANALYSIS 

W.  Hereman,  P.  P.  Banerjee  and  D.  Faker 

1 .  Introduction 

Over  the  last  five  years,  Painleve  analysis  of  ordinary  (ODEs)  and  partial 
differential  equations  (PDEs)  has  been  successfully  applied  to  get  an  idea 
whether  or  not  the  given  equation  or  system  -of  equations  might  be 
integrable.  In  essence  the  Painleve  test  verifies  if  the  solutions  of  the 
differential  equations  are  free  of  movable  critical  points;  in  other  words, 
the  only  movable  (i.e.,  dependent  on  the  initial  conditions)  singularities 
should  be  poles . 

With  reference  to  systems  of  nonlinear  coupled  ODEs,  the  Painleve  test 
provides  insight  into  the  possibility  of  constructing  invariants  (first 
integrals)  and  their  nature.  Direct  application  of  the  Painleve  test  at  the 
PDE  level  provides  an  algorithm  to  construct  the  infinite  set  of  conserved 
quantities.  Furthermore,  Painleve  analysis  serves  as  an  elegant  tool  to 
construct  the  Lax  pair,  and  thus  allows  to  linearize  and  subsequently  solve  the 
PDE  by  the  Inverse  Scattering  Transform  (1ST) .  Painleve  analysis  also  leads  to 
non-standard  auto-Backlund  transformations  by  suitable  truncation  of  the 
Laurent  expansion  of  the  solution.  For  ample  references  on  those  aspects  the 
reader  should  consult  the  papers  of  Weiss  (1983-1987) . 

In  this  paper  we  will  focus  on  yet  another  spin-off  of  the  Painleve 
analysis  of  PDEs.  More  precisely,  we  will  investigate  if  the  Painleve  test 
helps  to  construct  a  hierarchy  of  solitary  wave  solutions  to  the  Korteweg-de 
Vries  (KdV)  equation,  which  is  the  ubiquitous  equation  in  soliton  theory. 
Indeed,  Painleve  analysis  leads  in  a  natural  way  to  the  classical  (regular  and 
singular'  solitary  wave  solutions.  Unfortunately,  the  Painleve  approach  does 
not  permit  to  go  beyond  the  known  closed  form  traveling  wave  solutions 
to  the  KdV.  Although  the  associated  Painleve-Backlund  (PB)  equations  are 


invariant  under  two  distinct  classes  of  Schlesinger  transformations,  the 

application  of  the  latter  does  not  lead  to  new  results.  For  instance,  the 

Mobius  transformation  maps  solutions  of  the  PB  equations  to  solutions  to  the 

same  equations  but  with  different  values  of  the  parameters.  Within  the  Mobius 

group  we  only  find  closed  form  solutions  to  the  KdV  which  are  disguises  of  the 
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well-known  sech  and  cosech  solutions.  Application  of  the  Schlesinger 
transformations  to  rational  solutions  is  much  more  fruitful.  For  the  KdV  and 
many  other  equations  an  infinite  family  of  rational  solutions  can  be 
constructed  in  this  way.  This  topic,  which  is  extensively  discussed  and 
referenced  by  Weiss  (1983-1987)  is  outside  the  scope  of  this  paper. 

The  paper  is  organized  as  follows:  In  Section  2,  we  briefly  review  the 
basics  of  Painleve  analysis  and  derive  the  PB  equations  for  the  KdV,  which 
serves  as  our  paradigm  throughout  the  paper.  We  then  discuss  in  quite  some 
detail  the  properties  of  the  invariant  quantities  (such  as  the  Schwarzian 
derivative)  occuring  in  the  PB  equations.  The  most  important  result  is  that 
these  quantities  themselves  either  satisfy  the  KdV  or  the  modified  Korteweg- 
de  Vries  (MKdV)  equation. 

Next,  in  Section  3,  starting  from  the  PB  equations,  we  construct  the 
explicit  (closed  form)  traveling  wave  solutions  to  the  KdV.  The  validity  and 
also  the  limitations  of  this  technique  are  rigorously  proved  in  a  couple  of 
theorems.  The  method  is  illustrated  in  an  example. 

In  Section  4,  we  turn  to  a  study  of  the  various  connections  of  the 
Painleve  approach  with  other  well  established  methods.  We  start  with  Hirota's 
formalism,  which  leads  uo  N  solitary  wave  solutions  whereas  the  Painleve 
technique  apparently  does  not.  Second,  we  compare  with  the  Lax  method 
which  allows  us  to  linearize  the  KdV.  It  is  shown  how  the  Schrodinger  equation 
and  the  evolution  equation  for  the  eigenfunctions  follow  from  the  PB  equations 


by  the  balance  of  nonlinearity  and  dispersion  (or  dissipation)  in  the  given 
PDE.  For  particular  values  of  j,  which  are  called  the  resonances  and  are 
denoted  by  r,  it  will  be  possible  to  take  u^(x,t)  arbitrary  provided  the 
remaining  terms  in  the  recursion  relation  vanish  (compatibility  condition) . 

More  details  on  the  Painleve  test,  including  the  results  for  various  PDEs,  may 
be  found  in  papers  by  Weiss  (1983-1987),  Gibbon  et  al  (1985),  and 
Hereman  and  Van  den  Bulck  (1988)  .  The  latter  contains  a  MACSYMA  program  that 
symbolically  performs  all  the  required  steps  for  the  P-test  for  a  given  ODE  or 
PDE. 

Whether  or  not  the  Painleve  test  allows  us  to  draw  strong  conclusions 
about  the  integrability  (in  the  context  of  1ST)  of  a  given  nonlinear  PDE 
remains  debatable.  As  a  matter  of  fact  there  are  counterexamples  available 
(Gibbon  et  al  1985,  Clarkson  1985,  1986,  Fokas  1987)  disproving  the  necessity 
and  sufficiency  of  the  Painleve  property  to  guarantee  complete  integrability  of 
PDEs,  We  stress  that  integrability,  by  no  means,  refers  to  the  possibility  of 
constructing  particular  analytical  solutions  of  PDEs,  which  we  will  explore  in 
this  paper.  To  say  the  least,  Painleve  analysis  serves  as  a  handy  tool  to 
investigate  the  existence  of  solitary  wave  solutions  of  nonlinear  evolution  and 
wave  equations,  which  may  then  be  constructed  by  any  existing  direct  method 
(for  a  review  see  Kereman  et  al  1986) .  More  interesting  to  us  is  the 
derivation  of  auto-Backlund  transformations  via  a  truncation  of  (1)  and  the 
subsequent  construction  of  the  Lax  pair  (Lax  1968)  for  the  PDE.  A  truncation 
of  the  Laurent  series  is  possible  provided  the  coefficient  of  the  constant 
level  term  is  a  solution  to  the  original  equation.  Furthermore,  the  function  ^ 
must  satisfy  a  set  of  PB  equations.  These  equations  are  quite  easy  to  solve, 
specially  since  they  are  transformable  into  an  equivalent  set  which  is 
invariant  under  the  (Mobius)  group  of  homographic  (Schlesinger) 


transformations . 


as  well.  Finally,  since  a  non-standard  auto-Backlund  transformation  is 


obtained  by  truncation  of  the  Laurent  expansion,  a  comparison  with  the 
classical  Backlund  transformation  method  seemed  in  place. 

The  methods  presented  in  this  paper  carry  over  to  a  large  class  of 
completely  integrable  evolution  and  wave  equations.  As  a  matter  of  fact,  a 
study  of  the  Kuramoto-Shivashinsky  (KS)  equation  along  the  same  line  of 
thought  is  presented  in  a  paper  by  Conte  and  Musette  (1988a)  in  volume  1  of 
these  Proceedings. 


2 .  Painleve  Analysis 

Weiss  et  al  (1983)  proposed  a  direct  procedure  to  test  whether  or  not  a 
partial  differential  equation  is  of  the  Painleve  type  (P-type) ,  i.e.,  its 
solutions  admit  only  poles  as  movable  singularities.  The  test  requires  the 
substitution  of  a  Laurent  expansion  for  the  solution  of  the  PDE  in  terms  of  a 
new  dependent  variable  say 


u  (x,  t)  =  4)^  (x,  t) 


Uj  (x,t)<(c’  (x,t) , 


(1) 


and  subsequent  verification  if  the  recursion  relation  for  the  coefficients 
u^(x,t)  has  self-consistent  solutions  in  terms  of  the  time  and  spatial 
derivatives  of  (}>.  We  require  <|)(x,t)  to  be  analytic  in  some  neighborhood  of  the 
singular  manifold  (4i(x,t)=0),  which  itself  is  supposed  to  be  non¬ 
characteristic  (<Ji^  ^  0)  •  The  exponent  N  of  the  singularity  is  determined 


Consider  the  KdV  equation  in  u(x,t): 

u  +  ttuu  +  8u_  =  0,  (2) 

t  X  3x 

where  the  (rescalable)  real  coefficients  a  and  5  are  left  in  for  easy 
comparison  with  the  literature.  By  now  it  is  well-known  (Weiss  et  al  1983, 
Weiss  1983,  1984a)  that  (2)  passes  the  P-test  and  that  the  expansion  (1)  (with 
N  =  -2)  can  be  truncated  at  the  constant  level  term  u^,  to  give 


u (x, t ) 


125 

a 


X 

,2 


'^2x 


+  U2(x,t) 


125 

a 


(Kn  4')  2jj  +  ^2  ’ 


(3) 


provided  itself  is  a  solution  of  (2),  and  ([»  satisfies  two  PB  equations 
simultaneously.  An  easy  way  to  rederive  these  equations  is  as  follows. 
Subtract  the  KdV  in  u^  from  (2);  this  gives 


(u-U2)t  +  f  [(u-u^)^  +  2u2(u-U2)]^  +  ^<^"^2’3x  ° 


(4) 


Next,  substitute  u-u^  from  (3)  into  (4)  and  integrate  once  with  respect  to  x. 

-1  -2 

Setting  the  coefficients  in  (j)  and  (j>  equal  to  zero  results  in 


({I  +  au_4».,  +  5<t>.  =  0, 

^xt  2^2x  ^4x 


(5) 


Vt  “^2'^’x  ^  ^^‘^x''’3x  ■  “  °- 


(6) 


-3  -4 

The  coefficients  of  (]>  and  (|i  vanish  identically. 


Upon  elimination  of  u^  between  (5)  and  (6)  and  subsequent  integration 


with  respect  to  x,  one  obtains 
C  +  5s  =  X(t) 


(7) 


where  the  integration  constant  X  may  still  depend  on  t,  and  where  (Conte  and 


Musette  1988a) 


and  where 


/.  \  '**3x  3 

^2x 

-  4,  -  2 

X 

^  X  ' 

(8) 


(9) 


denotes  the  Schwarzian  derivative  (Hille  1916).  Apart  from  the  basic  ratios  C 
and  S,  which,  as  is  easy  to  show,  are  both  invariant  under  the  (Mobius)  group 
of  real  fractional  linear  transformations. 


o  I  '  ad-bc  *  0, 
c  +  d(p 


(10) 


there  is  a  third  ratio 

_  ‘‘’2X 


(11) 


which  is  not  invariant  under  (10) ,  but  nevertheless  plays  a  crucial  role,  as  we 
will  see  later  on.  Indeed,  from  (6),  u^  may  be  expressed  as 


u^  =  -  ^  (C  +  45s  +  35e^) . 


(12) 


Using  (8)  and  (9),  one  first  calculates 


1  2 

E  =  S  +  -  E  , 

X  2 


(13) 


E  =  C.,  +  EC  +  CS  +  ^  E^C  , 

t  2k  X  2 


(14) 


and  later  on  expresses  the  relevant  derivatives  of  u^  in  terms  of  C,  S,  their 
derivatives  and  E.  Straightforward  algebra  reveals 


(u^)  =  -  i  (C  +  4SS  -  eS^ES.,  -  66^ES 

2  t  a  t  t  2x  X 


+  66esc  +  36e  C) , 


(15) 


(16) 


(S  +  2ES^  +  E^)  , 

2  X  a  X 

*3  C  O 

(u^)^  =  -  —  (S*  +  6SS  +  8ES  +  lOE  S 

2  3x  a  3x  X 

+  5E^S  +  2ES^  +  3E^) ,  (17) 

X  2x 

where  we  have  used  (7)  to  eliminate  the  spatial  derivatives  of  C  in  favor 
of  those  of  S . 

Combination  of  (15),  (16)  and  (17)  according  the  KdV  equation  (2)  gives, 

after  simplification, 

C  +  55s  -  65SC  -  58CS  -  CC  +  48^  S.  +  28^SS  «  0  (18) 

t  X  X  X  3x  X 

or  equivalently  with  (7), 

X 

-  Xs  +  38ss  +  8s,  =  -  tI  .  (19) 

t  X  X  3x  4o 

Since  C  and  S  are  dependent  through  <{i,  the  compatibility  condition  (|>  ,  = 

u  I  3x 

d),  implies  the  constraint 
3x,  t 

S  -  C,  -  2SC  -  CS  =0.  (20) 

t  3x  XX 

Using  (7)  once  again,  one  thus  obtains 

S^  -  XS  +  35SS  +  6S_  =  0  (21) 

t  X  X  3x 

KdV  equation.  This  remarj^able  result  is  due  to  Weiss  (1984a,  1986b)  and  was 
rediscovered  by  Conte  and  Musette  (1988b),  who  brought  it  to  our  attention. 
Since  (20)  was  derived  independently  of  (5),  (6)  and  (7)  it  is  clear  from 

comparison  of  (19)  and  (21)  that  X^  =  0 ,  This  observation  allows  us  to 
conclude  that  for  constant  X,  the  requirement  that  U2  satisfies  the  KdV  is 
equivalent  to  a  compatibility  condition  in  terms  of  the  function  (|».  After  all. 


this  should  not  surprise  us  since  our  non-standard  derivation  of  the  Painleve- 


Backlund  equations  (5)  and  (6)  was  contingent  upon  u^  solving  (2) .  It  now 
suffices  to  solve  (7),  with  constant  X.  Even  this  can  be  slightly  simplified  by 
realizing  that  (7)  is  Galilean  invariant:  If  4>(x»t;X)  denotes  the  solution  for 
arbitrary  X,  then  (|>(x,t;0)  ••  4i  (x-Xt,  t;X)  .  Next/  (5)  needs  to  be  solved  for  u^ . 
Then,  substitution  of  (j)  and  u^  into  (3)  leads,  at  least  in  principle,  to  a  new 
solution  u  to  the  KdV.  At  this  point,  one  might  wonder  if  C  and  E  satify  any 
particular  well-known  equations.  From  (7)  it  follows  that 

C  +  2XS  -  3CC  +  6c,  =  0.  (22) 

t  X  X  3x 

The  equation  for  E  can  be  drived  as  follows.  From  (7)  we  calculate 

C.,  =  -  (EC)  (23) 

2x  t  X 

using  the  compatibility  condition  4>  =  <t)-  .  Since  C  can  be  expressed  in 

dm  f  ^ 

terms  of  S  by  (7),  and  S  relates  to  E  by  (13),  we  obtain  from  (23), 

E  -  Xe  -  ^  6e^e  +  6e,  =  0,  (24) 

t  X  2  X  3x 

which  is  the  modified  KdV  (mKdV) ,  a  result  also  due  to  Weiss  (1986b), 

3 .  Solitary  Wave  Solutions  Constructed  from  the  Painleve-Backlund  Equations 
In  our  search  for  closed  form  solutions  to  the  KdV  equations  we  adhere 
to  the  following  strategy: 

(i)  Starting  from  a  traveling  wave  solution  ^  in  real  exponential  form,  we 
calculate  the  dispersion  law  from  (7);  X  enters  a  free  constant 
parameter; 

(ii)  Next,  we  obtain  u^  from  (5) ; 

(iii)  Substituting  <t>  and  u^  into  (3),  we  find  u; 

(iv)  Subsequently  applying  the  Mobius  transformation  (10)  we  obtain  a 


new  4)  and  repeat  the  above  steps. 


(ad-bc)  -  _  a+b<t> 
d  c+dcj) 


Theorem  2 


Starting  from 


^  =  exp()cx  -  (iJt  +  6 


with  0Q  constant,  application  of  the  Mobius  transformation 
.  a+bd) 

=  ^  (29) 

will  not  lead  to  a  more  general  solution  to  the  KdV  than  the  one  obtained  from 
the  special  case  c  =  1,  d  =  0,  a  and  b  arbitrary. 

Proof 

Since  C  and  S  are  invariant  under  the  Mobius  group  (29)  one  can  evaluate  these 
for  <|»  as  in  (28)  .  With  (7)  we  have 


j  -  X  +  J  5k  , 


As  seen  from  (12),  only  E  needs  to  be  evaluated  to  find  u^ . 


Under  the  Mobius 


transformation 


E  =  )c 


e 

(c-de  ) 

e 

(c+de  ) 


which,  remar)cably  enough,  no  longer  depends  on  a  and  b.  From  (12),  obviously 


u  =  -  ^  (X  +  ^  [£n(c+de®)]  . 

2  a  2  a  2x 


Analogously,  from  (3) , 


u  =  -  ^  (X  +  I  8)t^)  +  (tn(a+be®)]  . 

a  2  a  2x 


This  procedure,  although  beautiful  in  theory,  is  of  limited  practical  use,  as 
will  be  shown  in  the  next  theorems  and  the  subsequent  example. 


Theorem  1 

Suppose  u  and  u^  are  solutions  to  the  KdV  equation  (2)  for  which  the 
difference  is  expressible  in  terms  of  the  second  logarithmic  derivative  of  a 


singular  manifold  (J),  which  is  lii  .ced  to  u^  by  the  PB  equations  (5)  and  (6). 
a  +  bd) 

Then  w  =  - “T  will  satisfy  the  PB  equations  (5)  and  (6)  provided  we  replace 

c  +  d(J) 

U^  by  u. 


Proof 

We  will  first  show  the  result  for  the  special  case 


Regarding  the  desired  result,  suppose  that  y  satisfies 


\u  +  av„w  + 
^xt  2^2x 


5y 


4x 


=  0. 


(25) 


(26) 


From  (25)  we  convert  all  the  needed  derivatives  of  V  into  derivatives  of  ((), 
and  we  substitute  into  (26)  .  Using  (5)  and  (6)  we  eliminate  the  t  derivatives 
of  (j),  elegantly  introducing  u^  into  the  expression.  Finally,  we  calculate  v^ 
as 


.  126 
2  2a 


1»2x«  -  ^ 


=  u. 


(27) 


Remark  that  (5)  and  (6)  are  invariant  for  scaling  of  4>  and  addition  of  a 
constant  to  it  (<^  ->  c  +  dA) ,  likewise  (26)  is  invariant  upon  replacing  y  by 
S  +  By.  Hence,  starting  from  c  +  d(^  we  apply  the  first  part  of  the  theorem  to 


y  *  1/  (c+d(^)  and  construct 


It  is  apparent  that  (33)  and  (34)  are  structurally  the  same.  From  Theorem 
1  it  follows  that  any  further  application  of  Mobius  transformations  would  not 
lead  to  anything  more  general  than  (34),  which  itself  follows  from  (29)  for  the 
siit^lifying  choice  c  =  1,  d  =  0.  Observe  that,  apart  from  a  constant,  u^  and  u 
in  (33)  and  (34),  have  the  same  "singular  part"  as  u  in  (3).  This  peculiar 
result  is  not  unique  for  the  KdV.  Even  evolution  equations  such  as  the 
Kuramoto  -  Sivashins)cy  (KS)  equation,  which  does  not  pass  the  P-test  in  the 
strict  sense  (due  to  the  presence  of  complex  conjugate  resonances)  exhibits  the 
same  behavior.  Conte  and  Musette  (1988a)  have  shown  similar  results  for  the  KS 
equation,  which  greatly  inspired  us  to  prove  theorems  1  and  2  for  the  famous 
KdV  case. 

Example 

Let  us  start  from  $  in  (28)  and  calculate  u  according  to  the  first  three 
steps  outlined  in  the  beginning  of  this  section.  The  dispersion  law  follows 
from  (30) : 

(0  =  -  (X  +  6)t^)  )c  .  (35) 


From  either  (5)  or  (6)  one  obtains 


“2  =  -  i  ^  I 


(36) 


For  the  present  choice  of  the  singular  part  in  (3)  vanishes,  so  that  u  «  u^ • 
We  now  carry  out  the  calculations  for 

0 

-  a  +  be  ,  (37) 

Q 

which  results  from  a  special  Mobius  transformation  on  e  ,  hence  X  and  0)  (see 
(35))  remain  the  same.  Following  the  notation  of  the  proof  of  Theorem  1, 


-  -  (X  +  I  8)c^) 
a  2 


u  *  u 


(38) 


Denoting  the  solution  corresponding  to  by  v,  we  obtain 


V  -  ^  (enV)2^  + 


126lt^abe®  1,,  ,3.2, 

-  -(X  +  2  Sk  I  . 

a(a+be  ) 


.  „  .  b 

This  solution  depends  only  on  the  free  parameters  A,  k  and  d^,  the  ratio  l^l 
being  absorbed  in  0^.  For  a=0  one  obtains  the  trivial  constant  solution  to  the 
KdV  equation.  For  ^  >  0  we  retrieve 


38k^  .2  k.  ^  ®0. 

V  =  -  sech  — [x  +  (A  +  —  6k  )t  +  — ] 

Ct  2  ^  jc 


i(X  *  f 


with  —  =  exp  0„  , 
a  0 

b  2  2 

For  —  <  0,  the  sech  in  (40)  is  replaced  by  -cosech  ,  the  singular  solution 

to  the  KdV  equation. 

Let  us  conclude  this  section  by  pointing  out  that  an  equivalent  result  can  be 
obtained  by  starting  from 


Y  =  A  +  B  tanh  -  , 


where  A  and  B  depend  in  a  suitable  way  on  a,b,c  and  d  in  the  Mobius  group 

a+b  c+d  ,  b 

(29),  where  only  the  parmeters  — and  are  relevant  (—  can  be 

b  d  d 

scaled  out) .  Calculations  along  similar  lines  as  before  produce  an  elegant, 
though  equivalent,  representation  of  the  regular  (and  singular)  traveling  wave 
solution  of  (2) : 


2  B  +  A  tanh  - 
6 

A+B  tanh  — 
^  2 


One  easily  verifies  that  (39)  and  (42)  are  the  same  for  A  “  — (a+b), 

1  2 
B  “  — (b-a)  .  For  a  graph  such  that  0  <  v  <  1  we  selected  a  -  -2\  -  36)c  for 

which  (42)  reduces  to 


V(;)  =  (l-A^)sech^C 
(l+A  tanh 


(43) 


where  A-B/A,  C“2  is  symmetric  with  respect  to  = 

-tanh  ^A.  Fig.  1  shows  v(^)  for  A  =  ^  .  For  any  A  >  1  the  solution  is 
2 

singular  (cosech  -  type) .  At  this  point  one  might  wonder  if  it  is  possible  to 
construct  N  solitary  wave  solutions  from  the  PB  equations.  The  answer  is  no, 
and  this  issue  will  be  addressed  in  the  next  section. 

It  is  however  possible  to  construct  an  infinite  class  of  rational 
solutions  to  the  KdV  equation  starting  from  (5)  and  (6),  and  by  exploiting 
invariance  properties  of  the  related  equation  (7)  other  than  the  Mobius  group. 
Remar)c  that  (7)  is  invariant  under  the  reciprocal  derivative 

transformation  4)  =  ^  (Weiss  1984a,  1986a, b,  1987) .  Details  on  the  construction 

of  rational  solutions,  which  is  out  of  the  scope  of  this  paper,  may  be  found  in 
papers  by  Weiss  (1984a,  1986b),  Gibbon  et  al  (1985)  and  references  therein. 


4 .  Connection  with  Other  Methods 
4 . 1  Hirota's  method 

Hirota  (see,  for  instance,  Ablowitz  and  Segur  1981,  Matsuno  1984) 
constructs  the  N  solitary  wave  solution  to  the  KdV  equation  (2)  in  a  quite 
ingeneous  way.  Hirota  substitutes 


a  2x 


(44) 


into  (2),  hence  reducing  it  to  a  quadratic  equation  in  f(x,t) 


ff  ^  -  f  +  8ff .  -  48f  f,  +  38f^  =  Gf^, 

xt  X  t  4x  X  3x  2x  ' 


(45) 


which  in  turn  can  be  written  in  bilinear  form  as 


(D  +  6  D^)  f»f  =  Gf^, 

X  t  X 


(46) 


using  the  bilinear  operator  (Newell  1985) 


(47) 


where  y  plays  the  role  of  an  auxiliary  variable.  Similar  formulae  hold  for  D 
and  products  (powers)  of  and  Setting  the  integration  constant  G  in 

(46)  equal  to  zero,  the  N  solitary  wave  solution  is  then  obtained  as  follows. 
Upon  substitution  of  a  formal  expansion  for  f,  say 


OO  *  /  '  \ 

f  =  1  (x,t) 


(48) 


into  (46)  and  equating  the  coefficients  of  the  boo)c)teeping  parameter  e  order 
by  order,  one  arrives  at  a  perturbation  scheme  for  the  subsequent  F^^l 
Starting  from  a  sum  of  N  real  exponentials 


,(1) 


N 


N 


=  ,r,  exp  (k.x  -  (O.t  +  0.)  =.Z,  e 
D=1  :  :  3  D=1 


0. 

3 


(49) 


the  expansion  (48)  can  be  broken  off  exactly  at  level  N  (i.e. 

=  ...  =0).  Calculation  of  the  RHS  of  (44)  gives  the  explicit  form  of  u(x,t) 

As  an  example,  the  two  solitary  wave  solution,  which  we  will  denote 
(2) 

by  u  (x,t),  IS  obtained  by  taking  N=2  in  (49),  which  will  satisfy 

3 

the  first  equation  in  th  •>  perturbation  scheme  provided  (0^  =  Sk^  (dispersion 

(2) 

law) .  After  some  tedious  algebra  involving  the  calculation  of  F  and 
subsequent  substitutions  of  the  results  in  (48)  and  (44),  one  gets 


2  0  2  9  ”2 
2  (k^  sech'^  ^2  cosech^  — ) 


0: 


(2)  36  ,.2 

u  =  —  (k„ 
a  2 


"i’ 


01  02  2 
(ki  tanh  -  k2  coth  — ) 


(50) 


.  65  /V.2  _  ^^*^2  ~  +  ^2  cosh  ®1  +  >^1  coshe2] 

~dr  '*^2  *^1'  - -  * 

Qi+Ot  ®1~®2  2 

[  (k2-ki)  cosh  ( — 2 — >  +  (k2+)ci)  co3h( — 2 — 

(51) 

Moloney  and  Hodnett  (1986)  have  shown  that  the  2  soliton  solution  can  be 
decomposed  into  two  amplitude  and  phase  modulated  solitary  waves.  For 
completeness  we  recall  their  results  here: 

u^^’=  ^  [kj  aO^)  sech^  I  (0^  +  HO^)) 

+  k^  A(e^)  sech^  |  (G^  +  0(0^))],  (52) 

with  amplitude  and  phase  modulations 

0'  20i 

1  +  Bje  ^  +  A22® 

A(0  )  =  - ^ ^ -  (53) 

(l+e  ^) (1+Ai2e  ^) 


H  (0  , )  =  en 
1 


1  +  Ai2e 


1  + 


with 


12 


^  -  ’'2 


ki  .  k^ 


and 


(54) 


(55) 


k^  ®2  ^  *^2  ®1  “  ^^*^1  “  ^2^^' 

The  ambiguity  in  selecting  can  be  resolved  by  requiring  that  the  total 

area  under  each  of  the  pulses  in  (52)  is  preserved  and  constant. 
Mathematically  this  implies  that 


-B, 


k^-k2 


Vki.k2 


-2 


(57) 


for  which  clearly  (56)  ia  satisfied.  Motivated  by  the  similarity  in  the 
starting  points  (44)  of  Hirota's  method  and  (3)  for  the  Painleve  approach,  one 
may  search  for  underlying  connections.  More  specifically,  one  may  wonder 
whether  or  not  it  is  possible  to  construct  multisoliton  solutions  from 
Painleve-Bac)clund  equations  (5)  and  (6) .  Recall  that  these  were  derived  by 
setting  the  terms  in  negative  powers  of  (}>  independently  equal  to  zero,  whereas 
for  (45)  there  was  no  such  restriction.  Hence,  multiplication  of  (5)  with  (|) 
and  subtraction  of  (6)  yields 


)  -  <b  (b  +  5(b<|) . 

^xt  ^x^t  ^  4x 


■  '  "’x'  ■ 


(58) 


We  can  now  see  that  Painleve-Backlund  equations  imply  (45)  for  G=0 
provided  u^  =  0  and  (|»  is  exponential  in  nature.  For  (j)  =  a  +  b  exp(kx  -  (i)t  +  0^) 
(58)  is  satisfied  since  (5)  and  (6)  are  with 


“2  ■  a  'k  -  ■ 


(59) 


For  u^  =  0  we  recover  the  dispersion  law  (i)  =  6k  used  in  Hirota's  method  and  f 
and  (J)  coincide.  In  Hirota's  formalism  a  two  soliton  solution  is  obtained  with 


f 


1  +  e 


+  e 


f  ^-^2 

I  kjtk^ 


) 


+  0. 


(60) 


which  surely  satisfies  (45)  for  G*0,  but  could  never  fulfill  (5)  and  (6) 
simultaneously  for  u^ •  Hence  it  appears  impossible  to  construct  multi-solitary 
wave  solutions  from  the  PB  equations  for  a  truncated  Laurent  expansion.  We 
strongly  believe  that  even  the  complete  Laurent  series  (see  Weiss  et  al  1983) 


would  fail  to  reconstruct  (50)  . 


As  a  final  remark,  it  is  possible  to  find  simple  solutions  of  (45)  for  G^O . 


For  instance,  recalling  that  two  functions  f  which  differ  by  the  exponential 
exp  (a(t)x  +  P(t)),  with  arbitrary  a(t)  and  P(t),  will  give  rise  to  the  same  u 
through  (44), 

f  =  cosh  exp  (-M0^) ,  (61) 

with 

1  2 

0  =  k(x  +  Xt  +  —  8k  t)  +  0Q  , 

and 

(X  +  I  8k^) 

M  -  - 

245k 

satisfies  (45)  provided 

G  =  -(2>.-5k^)  (2X  +  38k^)  .  (62) 

The  corresponding  solitary  wave  solution  is  in  (40),  from  which  f  in  (61) 
was  obtained  by  two  integrations  according  to  the  RHS  of  (44)  , 

4 . 2  Lax  method  and  Inverse  Scattering  Transform  (1ST) 

The  Lax  pair  for  the  KdV  equations  (Lax  1968)  can  be  derived  from  the  PB 
equations  (5),  (6)  and  (7)  in  an  easy  way.  First,  we  solve  (6)  for  4*^  and 

substitute  this  into  (7),  yielding 


35 


3x 


1 

X 


(au2  + 


X) 


(63) 


Secondly,  differentiation  of  (|i^  with  respect  to  x  allows  to  eliminate 

(b  between  (5)  and  (6)  resulting  in 
xt 


a(D^U2)x  +  35<d2<D^^  -  ^ 


0. 


(64) 


Next,  following  Weiss  et  al  (1983),  we  substitute 


X 


(65) 


into  (63),  to  get  the  Schrodinger  equation  (first  Lax  equation) 

65V2^  +  (au^  +  X)V  =•  0.  (66) 

/ 

Substitution  of  (65)  into  (5)  and  (64)  and  elimination  of  the  nonlinear  term 
between  the  resulting  equations  gives  the  second  Lax  equation 


V 

t 


+  au  V  +  ^(u  )  V  +  45v, 

2  X  2  2  X  3x 


0. 


(67) 


Defining  the  Lax  operators 


L 


65  3 
a 


+ 


PeR, 


(68) 


B  = 


45 


dx 


X  ®  /  3*  ^  3 

3  ^  2  <^^2  3l^  ^  ^ 


we  recast  (66)  and  (67)  into  a  compact  form  LV  = 


(69) 

AV  with  A  =  -  X/a, 


and  =  BV,  One  can  easily  verify  that  the  KdV  equation  (2)  in  u^  is 
equivalent  to 


L^  =  [B,L]  =  BL  -  LB.  (70) 

This  result  was  first  obtained  by  Weiss  et  al  (1983),  and  Weiss  (1983)  in 
two  different  ways,  and  the  same  procedures  for  obtaining  Lax  pairs  apply  for 
a  large  variety  of  ODEs,  systems  of  DDEs  and  PDEs.  For  the  latest  results  we 
refer  to  Newell  et  al  (1987).  From  (65)  the  connection  between  the  Painleve 
function  ^  and  the  eigenfunction  V  in  the  Schrodinger  equation  is  clear. 


4 . 3  The  Backlund  transformation  method 


Once  the  Lax  pair  is  available,  it  is  rather  straightforward  to  obtain 
the  Miura  transformation  and  even  the  traditional  Backlund  transformation. 
Indeed,  if  we  introduce  a  new  dependent  variable  r(x,t)  by  the  Cole-Hopf 
transformation  (Newell  1985) 

d 

r(x,t)  =  ~  =  ■^  V)  (71) 

V  ox 

or,  equivalently. 

V  =  exp  J  r(s,t)ds,  (72) 

then  the  Schrodinger  equation  may  be  replaced  by 

■^(r  +  r^)  +  pu-  =  A,  (73) 

ax  2 


which  is  known  as  a  Miura  transformation  between  the  functions  r  and  u^ . 

Since  u^  satisfies  the  KdV  it  is  reasonable  to  expect  that  r  will  satisfy  the 
modified  KdV  (mKdV) .  To  show  this  we  first  rewrite  the  second  Lax  equation 
entirely  in  terms  of  r  and  u^ : 


-r 

t 


=  126r^r  +  125r  ^ 

X  X 

+  46r  +  a(u  )  r 

3x  2  X 


*  2*“2*2x' 


(74) 


Subsequently,  we  eliminate  u^  between  (73)  and  (74) ,  yielding 

r  -  68r^r  +  r  +  Sr  -  0,  (75) 

t  X  p  X  3x 

which  indeed  is  the  mKdV.  The  traditional  Backlund  transformation  now  follows 
from  a  simple,  but  clever  argument.  Observe  that  if  r  is  a  solution  to  (75), 


(corresponding  to  u^)  then  -r  is  a  solution  as  well  (corresponding  to 

another  solution  to  the  KdV) , 

and 


(-r^  +  r^)  +  pu^  =  A. 


a 


(76) 


Subtracting  (76)  from  (73)  we  obtain 


its- 


(77) 


where  q  =  u  q  =  u  ,  q  and  q  being  the  potentials.  Substituting  the  last 

X  ^  f  X  ^ 

result  back  into  (73) ,  we  get  the  classical  Backlund  transformation. 


^  ,  2A  a  2 

*  “j’x  ■  r  ■ 


(781 


A  similar,  but  more  complicated  relation  for  (q  -  follows  from  (74), 

or  a  slightly  more  symmetric  form  (75)  .  Wahlquist  and  Estabrook  (1973)  have  shown 
that  this  evolution  equation  together  with  (78)  then  form  a  completely 
integrable  Phaffian  system.  Integrability  being  assured,  it  suffices  to  only 
consider  (78) .  Combining  (71)  and  (77)  we  also  have 


q  -  q  = 


125 

a 


3x 


(£n  V) 


or 


^2 


-  u. 


125 

a 


(€n  V) 


65 

a 


2x 


using  (65)  . 


(79) 


(80) 


Keep  in  mind  here  that  u^  and  u^  are  not  arbitrary  and  independent  solutions 
to  the  KdV.  To  be  precise,  u^  is  related  to  r  and  u^  to  -r.  From  (71)  it 
then  follows  that  u^  corresponds  to  V  whereas  u^  corresponds  to  1/V.  Going 
back  to  (65),  u^  corresponds  to  (|>  whereas  u^  corresponds  to  F  defined  as 


X 


(81) 


The  Backlund  transformation  (80),  which  involves  E  =  ^2-x.^^x.  defined  in 
(11),  suggests  that  (7),  which  only  involves  (j)  and  its  partial  derivatives,  is 
invariant  for  the  change  (j)  ->  T.  This  result,  which  was  first  observed  by 
Weiss  (1984a),  is  easy  to  verify  by  direct  calculation.  Using  (81)  in 
combination  with  (80)  gives 


^2  ■  '^2 


65 

a 


2x 


(82) 


Furthermore,  under  the  change  of  the  dependent  variable  in  (81), 


>  r  ' 

X  X 


(83) 


C*  =  7"  =  Cp  +  25 (Ep)  , 

<t>  r  r  X 


(84) 


Sx  =  (Ea)  -  |(E.)^  =  Sp  -  2(Sp)  , 

<j)  4)x24)  r  Tx 


(85) 


such  that  (7)  remains  invariant.  The  solution  corresponding  to  (|>,  denoted  by 

u.,  =  (u_).  will  be  transformed  into  u-  =  (u.)-  according  to  (80)  or  (82). 

2  2  <p  2  2  1 

Remar)c  that  in  eqs .  (83) -(85)  and  the  discussion  above  the  subscripts  <|)  and  T 
do  not  represent  partial  derivatives.  For  a  further  discussion  of  the 
Bac)^lund  transformation  (78)  we  refer  to  Wahlquist  and  E3tabroo)c  (1973)  . 


5 .  Conclusion 

Painleve  analysis  of  the  KdV  leads  to  a  non-standard  Baclclund 
transformation  involving  a  (singular)  manifold  (|»  which  has  to  satisfy 
simultaneously  two  partial  differential  equations  in  which  the  other  solution 
u^  to  the  KdV  occurs  as  a  coefficient.  Upon  elimination  of  the  latter 
(i.e.  u^),  one  obtains  a  partial  differential  in  ^  which  has  two  terms,  the  ratio 
and  the  Schwarzian  in  <),  that  are  invariant  under  the  Mobius  group  and  a 


reciprocal  derivative  transformation.  These  invariances  allow  us  to  construct  a 


hierarchy  of  new  solutions,  at  least  if  the  initial  (ji  is  rational.  For  <|>  of 
exponential  form  it  is  impossible  to  go  higher  than  one  rung  on  the  ladder: 
the  Mobius  group  remains  invariant  up  to  a  trivial  change  of  constants.  The 
meaning  of  the  function  <|>  is  illuminated  through  the  connections  of  the 
Painleve  analysis  with  Hirota's  perturbation  method,  the  Lax  method,  the 
Inverse  Scattering  Technique  and  the  standard  Backlund  tansformation  method. 

Acknowledgments 

The  authors  are  grateful  to  the  organizers  of  the  WASDA  III  Conference 
for  their  hospitality.  The  authors  would  also  like  to  acknowledge  the  support 
of  the  National  Science  Foundation  under  Grant  no.  ECS-8603643,  the  National 
Foundation  Presidential  Young  Investigator  Award  under  Grant  no.  MIP  8658865 
and  a  matching  grant  from  Xerox  Corporation,  a  grant  from  the  Rome  Air  Force 
Development  Center,  and  from  the  Air  Force  Office  of  Scientific  Research  under 
Grant  no.  85-NM-0263. 


References 


Ablowitz  M  J  and  Segur  H  1981  Solitons  and  the  Inverse  Scattering  Transform 
SIAM  Studies  in  Applied  Mathematics  4  (Philadelphia:  SIAM) 

Clarkson  P  A  1985  Phys .  Lett.  A  109A  205-8 

Clarkson  P  A  1986  Phys.  D  18D  209-10 

Conte  R  and  Musette  M  1908a  Proc.  Waves  and  Soliton  Days  Antwerp  III  Vol . I 
(Antwerp:  UIA)  3-19 

Conte  R  and  Musette  M  1988b  Private  Communication 
Fokas  A  S  1987  Stud.  Appl .  Math.  77  253-99 

Gibbon  J  D,  Radmore  P,  Tabor  M  and  Wood  D  1985  Stud.  Appl.  Math.  72  39-63. 

Hereman  W,  Banerjee  P  P,  Korpel  A,  Assanto  G,  Van  Immerzeele  A  and 
Meerpoel  A  1986  J.  Phys.  A:  Math.  Gen.  19  607-28 

Hereman  W  and  Van  den  Bulck  E  1988  Proc.  of  Workshop  on  Finite 

Dimensional  Integrable  Nonlinear  Dynamical  Systems  eds  Leach  P  G  L  and  Steeb  W- 
H  (Singapore:  World  Scientific)  117-29 

Hille  E  1976  Ordinary  Differential  Equations  in  the  Complex  Domain  (New  York: 
Wiley) 

Lax  P  D  1968  Comm.  Pure  Appl.  Math.  ^  467-90 

Matsuno  Y  1984  Bilinear  Transform  Method  (New  York:  Academic) 

Moloney  T  P  and  Hodnett  P  F  1986  J.  Phys.  A:  Math.  Gen.  11  L1129-35 

Newell  A  C  1985  Solitons  in  Mathematics  and  Physics  CBMS-NFS  Regional 
Conference  Series  in  Applied  Mathematics  £8  (Philadelphia : SIAM) 

Newell  A  C,  Tabor  M  and  Zeng  Y  B  1987  Physica  D  29  1-68. 

Wahlquist  H  D  and  Estabrook  F  B  1973  Phys.  Rev.  Lett.  21  1386-90 

Weiss  J  Tabor  M  and  Carnevale  G  1983  J.  Math.  Phys.  24  522-6 

Weiss  J  1983  J.  Math.  Phys.  ^  1405-13 

Weiss  J  1984a  J.  Math.  Phys.  ^  13-24 

Weiss  J  1904b  J.  Math.  Phys.  25  2226-35 

Weiss  J  1985a  J.  Math.  Phys.  26  258-69 


Weiss  J  1985b  J.  Math.  Phys .  26  2174-80 

Weiss  J  1986a  j.  Math.  Phys.  ^  1293-305 

Weiss  J  1986b  J.  Math.  Phys.  22  2647-56 

Weiss  J  1987  J.  Math.  Phys.  22  2025-39 


